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The resonances of the earth ionosphere cavity are considered directly as a cavity problem rather 
than as a limiting case of propagation in a quasi-waveguide with an arbitrary (spherically symmetric) 
conductivity profile. By using an approximation to a theoretically derived profile, values for the 
frequencies and Q factors of the lower resonances are computed which compare with experimental 
results as well as or better than the values derived from previous models. The elf eels ol small changes 
in the profile are also considered. It is shown that ionic conduction in the lower atmosphere should 
not be neglected and that the "knee" in the conductivity profile where the transition from ionic to 
electronic conductivity occurs has a significant, and previously unrecognized, effect, particularly on 
the Q factors of the resonant modes. 



1. Introduction 

The existence of resonances of the cavity between 
the earth and the lower edge of the ionosphere and 
their excitation by lightning flashes was predicted by 
Schumann [1952]. Experimental measurements have 
amply confirmed Schumann's prediction, and I he 
resonant frequencies and Q factors have been deter- 
mined [Balser and Wagner, I960], The "quality 
factor," Q, is defined as 2tt times the ratio of the stored 
energy to the energy loss per cycle; it is a measure 
of the sharpness of the resonance. Theoretical cal- 
culations have been performed under various simpli- 
fying assumptions [e.g., see Schumann, 1952; Wait, 
1964; Galejs, 1961a, b, and 1962; Chapman and Jones, 
1964; Jones, 1964]; many of these give reasonable 
overall agreement with experiment, but none is entirely 
satisfactory. The most successful models tend to be 
somewhat artificial and to show great sensitivity to 
artificially introduced parameters. Careful adjust- 
ment of parameters has been particularly necessary 
to reproduce the observed increase of Q with increasing 
frequency. 

It is conventional to approach the resonances as 
limiting cases of propagation in the quasi-waveguide 
between earth and ionosphere and to express the prop- 
erties of the boundaries by surface impedances. The 
models that give the best results consider a uniformly 
conducting earth and either an exponential ionospheric 
conductivity profile, sharply bounded at some (often 
frequency-dependent) critical height [e.g., see Galejs, 
1961a, b 1962], or an ionosphere constructed of several 
homogeneous layers [e.g., see Jones, 1964]. The 
effects of finite conductivity below the ionosphere 
may be treated by a later perturbation calculation. 

Treatment of the problem explicitly as a resonant 
cavity and use of a smooth conductivity-height profile 
to eliminate the nonphysical surface charge at the 
lower "edge" of the ionosphere should simplify the 
analysis and produce better results. This has been 



done by using a conductivity profile recently derived 
theoretically from fundamental considerations by Cole 
and Pierce [1965]. 

2. Differential Equations 

Maxwell's equations and the constitutive relations 
in mks units are 

VxH = J+f D D = Ke E 

at 

VxE + |-B = B = /x„H 



dt 



V-B = 0. 



J = crE 



Here H is the magnetic field, E the electric field, B 
the magnetic flux density, D the displacement, J the 
current density, and p the charge density. The per- 
mittivity and permeability of free space are denoted 
by €o and p, , respectively; K is the local dielectric 
constant and a the conductivity. 

The wave equation in H is derived in the usual way 
by assuming a universal exp (icot) time dependence, 
where co may be complex, taking the curl of the V X H 
equation, and eliminating B, E, and D. We find 



V 2 H + jUo/£€oO> 2 H — ijJL()(T(j)\\ 



+ 



V(cr +;Ke«&>) x(VxH) 



cr+iKe 



0. 



<)CO 



Now, letting fjL () e {) = ~ 1 
tion becomes 



, k 



and 17 = — , the equa- 



tor 



V 2 + Kk 2 -iki)- 



V(tj + iKk) x V x 



V 



iKk 



H = 0. 
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For vertical polarization (the conventional case), 
H • r = 0, where r is the radius vector from the 
center of the earth, and, since V • H = 0, we may 
expand H in spherical coordinates as 

lm 

where Xf(0, <p) = ===== r X V Yf(6, <p), and the 

1 V/(/+i) ' 

Yf(d, <p) are the ordinary spherical harmonics. 

Substitution of this expansion into the differential 
equation, with the assumption that 77 and K are func- 
tions of r only, gives the differential equation for the 
radial function fim{r) fi 



\dr 2 



/(/-hi) 



+ Kk 2 -ikrj- 



dr 



(r) + iKk) 



d 



7} + iKk dr 



/Kr) = 0, 



where "ra" has been suppressed because it does not 
enter the operator. The necessary boundary condi- 
tions are that the energy flow should be outward from 
the cavity, downward into the earth, upward into the 
ionosphere. 

If we further specialize to the case n = 7) e ~ a constant, 
for r < a, the boundary condition at r= a may be made 
explicit in simple form. 

For r < a, within the earth, 7] ~ 10 5 km -1 , r ~ 6 X 10 4 
km, and \k\ ~ 10 -4 km -1 at ELF, where \k\ is the mag- 
nitude of (complex) /r, so the leading terms of the 
differential equation are 



dr 2 ' 



ikr\ e 



f(r) = 0. 



The appropriate solution, corresponding to a downward 
wave that decays with time, is 

/(r)=^(exp[(^T7e) 1/2 r], 

where the square root is chosen to be in the first 
quadrant for k in the first quadrant. 
Just below the earth's surface, 



dr 



f(a-) = (ikn e ) ll2 f(a-). 



Assuming the earth-air transition to be sufficiently 
abrupt that the -j- (t) + iKk) term in the radial equa- 
tion predominates, we may write 

<L(4f\ 1 41- d. {r} + iKk)/iv + iKk) 



dr \dr) / dr dr 



so that 

df 7){a + ) + ik df 7){a + ) + ik df 

— (a+) = , r= , "7" \CL-) ~ , \Q>-)' 

dr 7) e +iK e k dr 7] e dr 

Choosing A so that/(a_)= 1, we have the initial condi- 
tions at r=a+ 



f(a + )=l 

f'(a + ) = (7 1 (a + mk)(ikl7 }( ,yl 2 . 

The boundary condition for large r may be expressed 
by requiring that the outward radial component of 
the Poynting vector, which may be written as 

Sr = (EXH*) r 

--[(rj-Im A)(Re/Re/' + Im/Im/') 

+ Re A(Re / lm /' - lm / Re /'), 

remain positive. Here * denotes complex conjuga- 
tion, Re and lm denote real and imaginary part, respec- 

df 



tively, and/' = 



dr 



In the frequency range, ~ 7 to 100 c/s, appropriate 
to the lower Schumann resonances, \k\ and 1/r are of 
comparable magnitude. At low altitudes, 7) < <\k\; 
at high altitudes, 77 > > \k\. For likely conductivity 
profiles, equality occurs at heights of 50 to 60 km. 
The differential equation is explicitly soluble in these 
limits, if we assume an exponential profile, 770 exp {rjL) 
for 7) at high altitudes. The solutions are 

7)<<\k\ ((P/dr 2 + k 2 )f=Q 



7]>>\k\ 



f=A exp (ikr) + B exp (— ikr) 
((Pldr^-ikT]^ exp (r/L)- \\L dldr)f=0 



/-exp (r/2L)[C//^(2L(-/A^ ) 1/2 exp (r/2L) 

+ DH?\2L(-ifa)o) 112 exp (r/2L))] 

where H^ and H {2) are the Hankel functions of order 
one of the first and second kind. If the square root 
is taken in the upper half plane, H^\H^) represents 
the solution that is regular (irregular) at infinity. 
Only the first is appropriate to this problem; the 
second represents incident energy. A "matching" 
of the asymptotic solutions to find the eigenvalue k, 
and thus resonant frequency and Q, would be highly 
inaccurate; because both the regular and irregular 
solutions are constant to first order at the lower end of 
their domain of applicability, the transition region 
where 7] ~ \k\ is of great importance. 

For any k a solution — which does not in general 
satisfy the boundary condition at large r— may be 
found by numerical integration of the radial differential 
equation, starting from the known boundary condition 
at r=a+. For k not an eigenvalue, S r will become 
negative above some height; an inward flux of energy 
is necessary to support the imposed variation of H 
with <p, 6, and t. By repeated integrations with suc- 
cessive refinements of the trial value of k y systemati- 
cally chosen to maximize the height at which S r 
becomes negative, we may approximate the eigenvalue 
ki to any desired degree of accuracy. Restated, the 
height h c at which S r becomes negative, considered 
as a function of complex /c, has a pole at the eigen- 
value ki. We are finding this pole by a simple ascent 
method. 
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Having found ki, the resonant frequency and Q 
may then he computed as 

/= c Re kl2ir 

Q = Re A/(2 ImA:). 

The result for i) > > \k\ strongly resembles one 
given by Wait [1962] for propagating waves in the 
flat-earth case. If the dropped terms k 2 — /(/ + 1)1 r 2 
are retained, with r approximated by the earth radius 
a, the order of the Hankel function is [1 — 4L 2 (A 2 

— /(/ + l)/a 2 )] 1/2 . Since the modified Bessel function 

77 

K is defined by K v (x) = —i p + 1 H ( £\ix), our solution dif- 
fers from Wait's only in that k 2 is replaced by k 2 

— /(/ + 1)1 a 2 , a familiar average correction for spherical 
geometry. 

3. Calculations 

The conductivity profile derived by Cole and Pierce 
leads to a profile for 7) that is reasonably well fitted by 
an equation of the form r)(h) = A exp (hl<x) + B exp 
(A//3), where h = r — a; the error is certainly of less 
significance than that due to the assumption of spheri- 
cal symmetry. The values of A, B, a, and /3 derived 
from the suggested approximation give 

7](h) = 7.5 X 10~ H exp (A/6.4) + 2.3 

X 10- 12 exp (A/3.0) km '; 

nearly as good a fit may be obtained with somewhat 
different parameters. The two terms represent ionic 
and electronic contributions to the conductivity. 

The value of r\ ( . should be between 3 X 10 3 and 
2 X 10 6 km -1 ; computed from the conductivities of 
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average earth and sea water, respectively. A reason- 
able average is probably about 10 5 km -1 . 

The radial equation was numerically integrated by 
using a high speed computer optimizing the complex 
variable k by successive approximation. This was 
done for the 77-profile above, referred to as I, and for 
two others: 

11 tj = 5.0X !()-" exp (A/6.4) 

+ 2.3 X 10- 12 exp (A/3.0) km" 1 



III T) = 5.0X10- 8 exp (A/6.4) 

+ 2.3X10- 



exp (A/2.7) krrr 



Profile II gives reduced conductivity in the ionic 
region and displaces the "knee" of the curve, while 
111 has a slightly steeper onset of electronic conduc- 
tion with an essentially unmoved knee. The value 
r\ ( >— 10 5 km -1 was used in most calculations; values 
of 10 4 and 10 H were also used for profile I, / = 3, to 
check sensitivity to this parameter. 

These profiles and the resulting mode frequencies 
and Q factors, together with those of several of the 
better previous models of various types, and experi- 
mental data are given in figure 1 and table 1. 

Note that the experimentally observed increase in 
Q with increasing frequency is reproduced in each 
case, I, II, and III, which has not been true of pre- 
vious models. Also of interest is the considerable 
effect, particularly on Q but also on frequency, of a 
small change in the conductivity in the ionic region 
(profile II). In previous models, this conductivity 
has been omitted or, at best, treated as a perturbation 
affecting Q but not frequency. 

The results are quite insensitive to earth conduc- 
tivity: with profile I, / = 3, a factor of 10 change in 
r) r changes / by 0.2 percent, Q by 0.5 percent, both 
quite insignificant. 
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Figure 1. Conductivity profiles. 
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Table 1. Frequencies and Q factors of the cavity modes 





/=1 


1 = 2 


/ = 3 


/ = 4 


1=5 


Profile 


Measured 
(1) 


/= 7.8 
# = 4.0 


14.1 
4.5 


20.3 
5.0 


26.3 
5.5 


32.5 
6.0 




Lossless 


/= 10.6 


18.3 


25.9 


33.5 


41.1 




Constant h' 

(2) 


/=8.3 


14.2 


20.0 


25.8 


31.7 
" 5.8 


•n=1.6x 10 10 exp (A/3.25) km ~ l h>h' = 50 km 
0<h<h' 


cr{h') = o)€o 

(2) 


/'= 8.0 
8 = 3.8 


13.8 


19.5 


25.4 


31.3 
"4.8 


r}= 1.6 x 10- 10 exp (A/3.25) knr'/i > h' 


"Average" 

(3) 


/'=8.05 
»C> = 6.75 
a c 4.5 


14.25 


20.5 


26.8 




tj=1.6X 10" 10 exp (A/3.25) km 1 //' = 45 km (Day) 


h 5.7 

'"•4.7 


1?= 8.1 x lO" 16 exp (A/2.25) km 'A' = 60 km 'Night) 


Two Layer 

(4) 


/= 7.90 
# = 4.7 


14.13 
6.0 


20.3 
6.9 


26.4 

7.5 


32.6 
8.2 


17=10 0<h<70 km 
1 X 10 4 km ' 70 < A < 100 km 
[4.9 X 10 5 km" 1 A > 100 


I 


/"=7.43 
0=4.1 




19.4 
4.7 




31.5 
5.1 


t? = 7.5X 10" exp (A/6.4) + 2.3 X 10" l? exp 
T7>= 10 5 km ' (A/3.0) km ' 














f±0.2% 
Q±0.5% 






Same, -n,, = 10 fi , 10 4 


II 


/=7.55 
(? = 4.5 




19.6 
5.3 




31.9 
5.6 


tj = 5.0x 10 « exp (A/6.4) + 2.3 X 10 '-' exp 
17,. = 10' km 1 (A/3.0) km" 1 








III 


7=7.71 
C> = 4.6 


13.9 

5.1 


20.0 
5.4 


26.2 
5.6 


32.4 
5.8 


tj = 5.0X10 * exp (/j/6.4) + 2.2Xl0 1: ' exp 
i7,= 10 r ' km ' (A/2.7) km ' 



(1) Balser and Wagner (I960 1. 

(2) Galejs |1961b|. 

(3) Galejs [1962J. 

(4) Jones [1964|. 
; ' At f= 10. 

" At /= 30. 

' With "cosmic ray" ionization. 17 exponential below A' 77 = 3.76 X 10~ 8 br 1 at A = 0, and continuous at A'. 



A typical modeform is shown in figure 2. Although 
drawn for the 1=1 mode with profile I, these profiles 
are all but indistinguishable from those for other 
modes and conductivity profiles. Energy and loss 
density profiles derived from this are given in figure 3. 
The energy densities are § /XoH • H* and \ e E • E*; 
the loss density is |E • J*|. Note that the loss density 
is significant only near the height where T7=|A:| (foi 
Im k < < Re k this is co = cr/eo). 

4. Importance of the Knee of the 
Conductivity Profile 

From the definition of Q as 2tt times the ratio of 
stored energy to energy loss per cycle, we have 



The spherical harmonics are normalized to 

f X* -X!<ffl= f Y*Yi(Kl=l, 
so that (neglecting field penetration into the earth) 

co 1(1+1) f X |/| 2 /fo 2 +|*| 2 )dA + a 2 f °° \ft 2 dh 
n Jo J o 



2c 



Ki+i) r i/i 2 /(^+i^i 2 )^ 

oJ 



<?= 



a)[e„JE • EW+^o/H • H*dV\ 



2/E • PdV 
In the cavity, r » a, the earth radius, and 

H «^X,=//a(/(/- I))" 1 ' 2 ^ Y,(6, <p)e e 
E r - -Jla^+ie^)-' {—^ ctn dj^j {l{l+\)Y^Yi{0, tp) 



Qi 



1 

itL[ 



, (Re kid) 2 
/? '" + 1(1+1) h < 



cfi f 



ar T) + ik 



d(i+i)yime,<p) 



E ~ E r e r at low altitudes 

nf 



Jr = CrE r 



a 2 (rj + ik) 



W+Dmie,?). 



Now if we approximate r) by 770 exp (r/L) near h , where 
7}= \k\ , and observe that "|/| ~ 1 for < h < A*, and 
/~ elsewhere, we have upon performing the integrals 



where Re A ~ \k\ has been used. 

Note in figure 1 that the points where \k\=r) for 
1=1 to 5 (\k\ ~ 2 to 7 X 10- 4 km- 1 ) lie on the knee of 
the curve. This has two effects: the effective value of 
L decreases as / increases, and the increase in h as / 
increases is greater than it would be for a single ex- 
ponential fitting the upper portion of the conductivity 
profile. Both tend to produce an increase in Q with 
increasing /. Therefore neglect of ionic conductivity, 
which is responsible for the knee, is not justified. 
Treating it as a perturbation to compute corrections 
to the Qs, as has been done before [Galejs, 1962], is an 
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FIGURE 2. Approximate radial dependence of fields. 



improvement, but the changes in the resonant fre- 
quencies (clearly not negligible: compare cases I and 
II) are all but impossible to calculate. 

5. Conclusion 

The values of the resonant frequencies and of the 
Qs of the cavity resonances found by the method of 
this paper compare with experimental values as well 
as or better than those of previous models. Although 
further variation of parameters would almost certainly 
produce more accurate results, the worth of such cal- 
culations is questionable. First, the experimental 
values, particularly those for Qs, are not that well 
established, and the resonant frequencies exhibit 
considerable temporal variation. Second, the assump- 
tion of spherical symmetry is clearly invalid because 
the conductivity profiles change temporally (diurnally, 
seasonally, etc.) and with geographic location; these 
influences may be expected to have an effect that 
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FIGURE 3. Approximate radial dependence of energy and loss 
densities. 



should not be neglected. Also, anisotropy of the 
ionosphere due to the earth's magnetic field will have 
some effects. 
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